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Leo IV & V - A possible dwarf galaxy pair ? 
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ABSTRACT 

The last few years have seen the discovery of many faint and ultra-faint dwarf spheroidal galaxies around the Milky Way. Among 
these is a pair of satellites called Leo IV and Leo V. This pair is found at large distances from the Milky Way (154 and 175 kpc 
respectively). The rather small difference in radial distance, and the fact that they also show a close projected distance on the sky, has 
led to the idea that we might be seeing a new pair of bound galaxies - like the Magellanic Clouds. In this paper we investigate this 
speculation by means of a simple integration code (confirming the results with full N-body simulations). As the luminous mass of 
both faint dwarfs is far too low to allow them to be bound, we simulate the pair assuming extended dark matter haloes. Our results 
show that the minimum dark matter mass required for the pair to be bound is rather high - ranging from 1.6x10'" Mq to 5.4 x 10'" M© 
(within the virial radii). Computing the mass of dark matter within a commonly adopted radius of 300 pc shows that our models are 
well within the predicted range of dark matter content for satellites so faint. We therefore conclude that it could be possible that the 
two galaxies constitute a bound pair. 

Key words, galaxies: dwarf - galaxies: binaries - galaxies: halos - galaxies: individual (Leo IV, Leo V) - methods: numerical 



1. Introduction 

The last decade has seen the discovery of many new faint 
dwarf spheroidal (dSph) galaxies of the Milky Way (MW) 
(e.g. W illman et aL, 2005; Belokurov et al., 2006; Zucke r et al. . 
I2006al; iBelokurov et al.1 120071; IWalsh et al.L |2007, and many 
more). Many of these dwarfs are less luminous than a globu- 
lar cluster (or even an open cluster) and exhi bit high velocity 
dispersions (given their luminous mass) (e.g. ISimon & Gehal 
l2007; Koch et al., 2009; Geha et al., 2009). Should these ob- 
jects be in virial equilibiium, they are the most dark matter 
(DM) dominated objects known in the universe. They would ex- 
hibit mass-to-light (M/L) ratios of rnore than a thousa i id (e.g . 
ISimon & GehaL I2007t iFellhauer et al.L l2008t iGeha et al.L 12009^ 

ACDM simulations (e.g. Millennium II simulation of 

Bovla n-Kolchin et al., 2009), (Via Lactea INCITE simulation of 
Kuhlen et al., 2008) predict that a galaxy like our MW should be 



surrounded by hundreds, if not thousands of small DM haloes 
which could host a dwarf galaxy. The discrepancy between the 
known number of MW satellites and these pre dicted values 
is known as the 'm issing satellite problem' (e.g. iKlvpin et al.L 
Moore et al. 1 11999'). The discovery of new faint dwarfs 
in the SDSS catalogue doubled the number of known satellites. 
Extrapolating to areas of the s ky an d distances not cove red by 
the survey (e.g. lKoposov et all 120081; [Maccio et aUl2009ll2010l) 
may suggest that the missing satellite problem is now solved. 

Amongst these ultra- faint galaxies are dwarfs which 
are extremely metal-poor (iKirbv et al.L l2008l and references 
therein), som e show complex star formation histoiies (Leo T 
Ide Jo ng et al. . 2008) and many o f them sho w unusual mor- 
phologies (Coleman et al., 2003 ISandetaLl 12009) Others 



seem to show signs of tidal disniption (IZucker et al.L l2006bl; 



IFellhauer et all 120071: iMunoz et aLlhoOSh . 



Signs of tidal disruption and the fact that all known dwarfs 
seem to be aligned in a disc-like structure around the MW (e.g. 
iMetz et all l2008l 120091) has given rise to alternative explana- 



Table 1. Observation al properties of Leo I V and Leo V. The data 
is mainly taken from lde Jong et al.] (l2010l) . 



Galaxy 



Leo IV 



Leo V 



RA (J2000) 


IP 32'" 58?6 


± 1?6 


11" 31™ 


8?4 ± 1^6 


Dec (J2000) 


00° 33' 6" 


±54" 


02° 12' 57" ± 12" 


Dist [kpc] 


154 ±5 




175 ±9 


vgsr [kms"'] 




10.1 




60.8 


Lv [Lo] 


1.8 ±0.8 


X lO'* 


1.0 ± 


0.8 X lO"* 


Mv 


-5.8 


±0.4 




-5.2 ±0.4 


'"h [pc] 




128 




133 


CTioj [kms-'] 


3.3 


± 1.7 




2.4-3.7 



tions for the existence of those galaxies. These theories imply 
that most, if not all, of the MW dwarfs are in fact no more than 
dis rupting star clusters or ti dal dwaif galaxies hosting no DM at 
all (iPawlowski et al.Ll201lh . 

In this paper we foc us on two of these ne w ultra-faint 
dwarfs, namely Leo IV (IBelokurov et al.L l2007l) and Leo V 
(IBelokurov et al.L |2008). The ir properties h ave been stud ied by 
many authors: for Leo IV: (iMoretti et al.l 12009; Siinon et al.L 
120 lOL ISand et al.L 1201 Ol) and for Leo V: (IWalker et al.L l2009allbl) . 
We summarize a selection of the observational data in Tab. [T] 
The pair of galaxies is found at a rather large distance from 
the MW (154 and 175 kpc respectively). The two galaxies are 
not only very close to each other in rad ial distance (22 kpc) 
(IMoretti et all l2009t IdeJong et al.L l2010l) but also in projected 
distance on the sky. Their radial velocity differs only by about 
50 kms In the discovery paper of Leo V the authors spec- 
ulate that the two dwarfs could be a bound pair similar to the 
Magellanic Clouds. In particular, the smalle r dwarf , Leo V , 
exhibits deformed elongated contours (Walk er et al.L l2009al) . 
which could be signs of tidal interaction. Furthermore, th ere is a 
tentative stellar bridge between the two satellites (de Jo ng et al.L 
2010), which now i s more hkely a f o reground stream of the 
Virgo overdensity d.Tin et alJ. l2012h . Ide Jong et all (l2010l) ar- 
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Table 2. Initial positions and velocities in the adopted Cartesian 
coordinate system. The final set of velocities assumes the same 
amount of relative tangential velocity as given by the relative 
radial velocity. 





Leo IV 


LeoV 


A 


X [kpc] 


14.713 


20.937 


6.224 


Y [kpc] 


84.728 


90.435 


5.707 


Z [kpc] 


128.442 


149.262 


20.840 


Vy, [kms-'l 


0.961 


7.242 


6.281 


V'y [kms-'] 


5.537 


31.282 


25.745 


[kms-'] 


8.392 


51.630 


43.238 


V^' [kms-'l 


-4.409 


39.276 


43.686 


[kms-i] 


12.945 


-14.784 


-27.731 


[kms-'] 


4.118 


75.047 


70.928 



gued that, to form a bound object, the twin system would need 
a lot of DM - much more than is seen in similar faint satellites. 
Nevertheless, the authors claimed that it is highly unlikely that 
the two satellites are a simple by-chance alignment. They also 
rule out the possibility that the two faint dwarfs are not galaxies 
at all but simple density enhancements of a stellar stream by or- 
bital arguments. They conclude that they might be a 'tumbling 
pair' of galaxies that have fallen into the MW together 

With our paper we want to investigate the hypothesis of a 
bound pair further Using a simple two-body integration method 
we investigate a large part of the possible parameter space, 
searching for the minimal total mass the system needs, in or- 
der to form a tightly bound pair We assume that both galaxies 
have their own DM halo, orbiting each other We describe the 
setup of our simulations in the next section. We then report and 
verify (using full N-body simulations) our results in Sect. [3] and 
finally discuss our findings in Sect.lH 

2. Setup 

2.1. Parameter Space 

We use the p ositions and distances, as rep orted in the discov- 
ery papers of iBelokurov et all (l2007l l2008l) . This enables us to 
transform the positions into a Cartesian coordinate system. The 
coordinates for Leo IV and Leo V are shown in the first three 
lines in Tab.|2] 

Current observations only provide radial velocities for the 
pair. Nothing is known so far about their proper motions. We 
therefore adopt two velocity cases in our simulations. First, we 
investigate the case that the radial velocity is the only velocity 
component the galaxies have (line 4-6 in Tab. ^ and second, 
we adopt an outward tangential velocity of the same magnitude 
as the radial velocity (line 7-9 in Tab. |2]i. Here 'outward tan- 
gential' means that the velocity vector is perpendicular to the 
radial velocity and points away from the other dwarf. The latter 
choice could be regarded as the worst case scenario. With these 
assumptions we restrict the possible parameter space of initial 
conditions significantly, but we are able to deduce how strong 
the influence of a tentative tangential velocity is on our results. 

To restrict our parameter space further we adopt just two 
cases for the mass-ratio between the two satellites. The observed 
absolute magnitudes give a luminosity ratio of 1.8. Adopting 
a mass-follows-light scenario for one of our cases, we use the 
same ratio for the two DM haloes: 



i-LeoIV 1-8 ■ 10 



iLeoV 1 ■ 104 



1. 



LeoIV 
LeoV 



Ml 



DM 



IWalker et al.l (l2009bh suggest that the faint and ultra- faint dwarfs 
reside in similar DM haloes of a certain minimum mass. 
Therefore, we investigate equal mass DM haloes as the other 
case. 

For the four cases described above we search for solutions 
adopting halo concentrations of c = 5, 10 an d 20, as those are 
values typically adopted for dwarf galaxies (Loka s & Mamoni 
I2OOI). This gives a total of 12 different solutions to the problem. 

The haloes are described as NFW-profiles with c - r^i^jr^^ii^, 
'"vii = '"200 being the virial radius in which the density is 200 
times the critical density of the universe, using a standard value 
of the Hubble constant of Hi) - 70 km s ' 



Mpc" 



2.2. Method 



Ide Jong et all ( 1201 01) investigated the minimum mass for the two 
DM haloes of the satellites, by assuming they were point masses. 
We use a simple two-body integration programme modelling the 
system in the following way: 

- Both satelhtes are r ep resent ed by analytical, rigid 
iNavarro. Frenk & Whitel (1 19971 NFW) potentials. The 
force on the centre of mass of one halo is computed using 
the exact force according to its position with respect to the 
potential of the other halo. 

- To be a tightly bound pair we adopt a rigid distance criterion 
which requires that neither centre of the two haloes leaves 
the halo of the other dwarf, i.e. their separation is always 
smaller than the virial radius of the (smaller) halo. 

- For each case we choose the total mass of the system and 
set up the two haloes according to their mass-ratio and 
concentration. Then we run the two-body code to determine 
if our distance criterion is fulfilled. If the maximum distance 
is larger or smaller, we alter the the total mass respectively 
and use the code again. We iterate this process until we find 
a maximum separation equal to our distance criterion. 

The reason why we choose a distance criterion instead of 
computing the escape velocity (i.e. the velocity the two dwarfs 
need to separate from each other to infinity) is that if we were 
to adopt such a criterion, we would also include bound cases in 
which the maximum separation between the two dwarfs easily 
exceeds their distance to the MW. As this would not make sense, 
we exclude these solutions by imposing a very rigid distance 
criterion. 

Furthermore, an estimation of the tidal radii of Leo IV and 
V gives values that are similar to our distance criterion. Using 



'"tidal 



'Wdwarf 

3Mmw(£>) 



1/3 



D 



(2) 



LeoV 



(1) 



with mdwarf ~ 4 X 10'" Mq, Mmw(£>) ~ 10'^ Mq and D * 
165 kpc, we get 40 kpc as tidal radius. This is slightly lower 
than the distance criterion used, but is also a rather rough esti- 
mate. However, it shows that the distance criterion is a sensible 
way to restrict the solutions. 

Hence, the choice of our distance criterion allows us to treat 
the galaxy pair as isolated, i.e. to neglect the potential of the 
MW. 

The simulations are always computed forward in time, start- 
ing from our current view of the dwarfs. We thus ascertain the 
next maximum separation to assess if the two dwarfs are bound 
to each other now. As we do not know where they came from, 
nor the details of their orbit around the MW, we cannot predict 
their future fate using these models. 
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Table 3. This table shows the minimum bound mass for each of our cases. The first column gives the number of the case, the second 
is the adopted concentration of the haloes. Then we give the mass of the DM halo and its virial radius for Leo IV and Leo V. The 
next column gives the total mass in DM of the whole system, the next column shows the 'ratio' by which the maximum distance 
differs between the full N-body simulation and the two-body code. The last column is a short explanation for the cases: rad. vel. = 
only radial velocity, rad. & tan. vel. = radial and tangential velocity adopted; mass ratio L8 = the two haloes have a fixed mass ratio 
of 1.8; equal mass = the two haloes have the same mass 



Case 


c 


Mdm.LcoIV 


vir,LeoIV 


Mdm.LcoV 


vir,LeoV 




ratio 


Scenario 






[Mo] 


[kpc] 


[Mo] 


[kpc] 


[Mo] 






Oa 


oo 


- 


- 






4.18 X 10" 


— 


Point masses (a) 


Ob 


oo 










1.47 X 10'" 




Point masses (b) 


1 


5 


1.34 X 10"' 


49.00 


7.45 X W 


40.29 


2.09 X 10'" 


0.965 


rad. vel. mass ratio 1.8 


2 


10 


1.27 X 10'" 


48.11 


7.05 X lO"* 


39.55 


1.98 X 10'" 


0.953 


rad. vel. mass ratio 1.8 


3 


20 


1.22 X 10'" 


47.42 


6.75 X W 


39.98 


1.89 X 10'" 


0.935 


rad. vel. mass ratio 1.8 


la 


5 


9.05 X 10'' 


42.99 


9.05 X 10'' 


42.99 


1.81 X 10'" 


1.179 


rad. vel. equal mass 


2a 


10 


8.55 X 10' 


42.18 


8.55 X 10'' 


42.18 


1.71 X 10'" 


1.194 


rad. vel. equal mass 


3a 


20 


8.30 X 10' 


41.76 


8.30 X 10'' 


41.76 


1.66 X 10'" 


1.204 


rad. vel. equal mass 


4 


5 


3.47 X 10'" 


67.25 


1.93 X 10'" 


55.28 


5.39 X 10'" 


0.993 


rad.& tang. vel. mass ratio 1.8 


5 


10 


3.11 X 10'" 


64.83 


1.73 X 10'" 


53.30 


4.83 X 10'" 


0.956 


rad.& tang. vel. mass ratio 1.8 


6 


20 


2.84 X 10'" 


62.90 


1.58 X 10'" 


51.70 


4.42 X 10'" 


0.968 


rad.& tang. vel. mass ratio 1.8 


4a 


5 


2.40 X 10'" 


59.50 


2.40 X 10'" 


59.50 


4.80 X 10'" 


1.450 


rad.& tang. vel. equal mass 


5a 


10 


2.20 X 10'" 


57.80 


2.20 X 10'" 


57.80 


4.40 X 10'" 


1.372 


rad.& tang. vel. equal mass 


6a 


20 


2.10 X 10'" 


56.91 


2.10 X 10'" 


56.91 


4.20 X 10'" 


1.284 


rad.& tang. vel. equal mass 



2.3. Full N-body simulations 

To ensure that the results are reasonable, we perform full N- 
body simulations as a check of each of the 12 solutions ob- 
tained with our simple code. W e use the particle-mesh code 
SUPERBOX (fellhauer et all 120001 It is fast and enables sim- 
ulations of galaxies on normal desktop computers. It has two 
levels of higher resolution grids, which stay focused on the sim- 
ulated objects, providing high resolution only in the areas where 
it is needed. 

Each object (halo) is modeled using 1,000,000 particles. We 
use NFW distributions for the haloes according to the results 
obtained with the two-body code. The haloes extend all the way 
to their virial radius. 

The resolution of the grids is such that we try to keep about 
15 cell-lengths per scale length, rjc, of the haloes. A particle- 
mesh code has no softening-leng th hke a Tree-code but previ- 
ous studies (Spin nato et al.L l2003 h showed that the length of one 
grid-cell is approximately the equivalent of the softening-length. 
Furthermore the particles in a particle-mesh code are not stars or 
in our case DM-particles, they rather represent tracer-particles 
of the phase-space of the simulated object. Densities are de- 
rived on a grid, and then a smoothed potential is calculated. 
The number of particles is chosen according to the adopted grid- 
resolution to ensure smooth density distributions. A detailed dis- 
c ussion about the partic le-mesh code SUPERB OX can be found 
in lFellhauer et al.l OOOOl) . 

To clarify we state once more, that the N-body simulations 
are only used to verify the results of the simple method. This 
means we check the next maximum distance of the two galax- 
ies. The simulations do not represent full-scale simulations of 
the past, present and future of the two Leo galaxies. Such simu- 
lations are more demanding and are beyond the simple scope of 
this paper. 

3. Results 

3. 1. Point Mass Case 

We first recalculate the minimal bound mass assuming both 
galaxies are point masses, following the methodology by 



6x10'° 
4xlO'° 



_ 1 ' 




1 1 1 1 




















o 








• 

1 < 




, 1 , , 


, , 1 



0.1 0.2 

1/c 

Fig. 1. Minimum DM mass of the total system Mtot versus con- 
centration of the haloes. We plot 1/c in favour of c to include 
the point-mass results at 1/c = (cases Oa,b; plotted as open 
and filled circle - blue online). Solid lines are the results us- 
ing the mass-ratio of 1.8. Triangles are cases with radial veloc- 
ity only and tri-stars with additional tangential velocity. Dashed 
lines show the results of equal mass haloes. Squares are radial 
velocity only cases and crosses have additional tangential veloc- 
ity. 



Ide Jong et all (12010 ^. We assume that the relative velocity be- 
tween the two dwarfs is equal to their escape velocity (i.e. the 
velocity required for the haloes to separate to infinity). Because 
point masses do not have any characteristic radius we cannot ap- 
ply any meaningful distance criterion here. 

Using only the observed radial velocities (cas e Oa) we con- 
firm the total DM mass of the system obtained bv lde Jong et al.l 



3 



M. Blana etal.:Leo IV & V 



Table 4. In this table we show the mass-range of the haloes of 
cases l(a)-3(a) and 4(a)-6(a) (which encompass the range of 
our results) within a radius of Leo IV of 97 pc, f or which an ob- 
serva tionally based estimate of the mass exists (ISimon & GehaL 
120071) 



Cases Km ('■"P^ = 9Vpc) M/L(r,,,) 

[Mo] [Mo/Lq] 



Simon & Geha 2007 


(1.4 ± 1.5) X 10" 


151 


1-3 


6.76 X 10^-4.60 X 10* 


37 - 255 


la -3a 


5.92 X 10^- 4.03 X 10* 


33 - 224 


4-6 


9.31 X 10' -6.19 X 10* 


52 - 344 


4a -6a 


8.23 X 10'- 5.58 X 10* 


46-310 



(I2010h . Our result is a factor of two lower, however, those au- 
thors used an approximation obtained from energy arguments, 
while we perform the full escape velocity calculation. Given the 
wide range in possible results, as shown later in the text, and tak- 
ing the observational uncertainties into account, a difference of 
a factor of two is still a very good match. 

Furthermore, we calculate the minimum bound mass assum- 
ing the two satellites also have tangential velocities (according to 
Tab.|2]i. In this case (case Ob) the total mass required to keep the 
galaxies bound is 1 .47 x 10'" Mq. That is, three times more DM 
is required to keep them bound. While this may seem a large 
mass, we refer the reader to the following sections to put this 
result into context. 

The results are shown in Tab. |3]in the first two lines and in 
Fig. [Hat l/c = 0. 





1 1 1 1 1 1 1 1 1 




' ' ' 
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LgoV 








, iw'' 


T 








Leo IV 










T 










UMa 












A 




- 






LeoV 


- 








T 








Draco Leoiv'' 
















- 




UMi*" ™« 




- 


- 




^Carina 




- 
















Sextans 




















Leoll 








Leol And II 










Fornax " * 










Scu Iptor 










1 , 1 1 , 1 1 1 









-14 -la -10 -B -6 -4 

Fig. 2. Mass-to-light ratios within a radius of 300 pc. The cir- 
cles are the d Sph galaxies of the Local Group as reported by 
I Wilkinson et al. (2006). For Leo IV, the filled upside-down tri- 
angles are cases la and 3 (concentrations c = 5 and c = 20 
respectively) with radial velocities only. The open triangles are 
cases 4a and 6 (c - 5 and c = 20) with radial and tangential ve- 
locities. For Leo V we plot the values of cases 1 and 3a as well as 
4 and 6a, respectively. We use these specific cases as they span 
the whole range of our results. 



3.2. Two-body Integrator 

We now have a look at our results, obtained by the two-body 
integrator we use. As explained in Sect. |2] the two haloes are 
rigid, analytical potentials acting on the centre of mass of the 
other galaxy. Additionally, we now introduce a very strict dis- 
tance criterion of the form that neither halo centre should leave 
the other halo (i.e. separations larger than the (smaller) virial ra- 
dius). This way we make sure that we are really dealing with a 
tightly bound pair 

The results are two-fold: Of course we see an immediate 
large increase in the required minimum mass (compared to the 
point-mass cases) just by introducing the rigid distance criterion. 
We plot the total mass in DM against 1/c in Fig.[T] to include the 
point mass cases (c - oo). We see that the bound mass is larger 
for lower values of the concentration. This can be easily under- 
stood as with higher concentrations we have more of the total 
mass of the halo concentrated towards the centre and therefore 
the gravitational pull on the other dwarf is larger 

Secondly, we also find that including the additional tangen- 
tial velocity roughly triples the required mass. While the cases 
with radial velocities (cases 1-3, la-3a) require masses in the 
range of about 1.6-2.1 x 10'" Mq, the additional tangential ve- 
locity increases the necessary masses up to much larger values 
of 4.2-5.4 X 10'° Mo (cases 4-6, 4a-6a). 

As a little side-remark, we see that we need slightly less 
massive haloes in the equal halo mass cases than if we adopt 
a mass-ratio of 1 .8 between the two haloes. As these differences 
are small compared with the differences of the unknown concen- 
tration, and even more so with the unknown tangential velocity, 
we can easily neglect them and assume that distributing the mass 
differently between the two haloes has no strong effect on our re- 
sults. 



Table 5. Mass-to-Light ratios within a radius of 300 pc of our 
simulations. We adopt V-band magnitudes of -5.8 for Leo IV 
and -5.2 for Leo V. 







LEO IV 


LEOV 


Case 


c 


Logio(Mo/Lo) 


Logio(Mo/Lo) 


1 


5 


2.544 


2.710 


2 


10 


2.929 


3.092 


3 


20 


3.343 


3.499 


la 


5 


2.484 


2.740 


2a 


10 


2.867 


3.122 


3a 


20 


3.279 


3.534 


4 


5 


2.686 


2.853 


5 


10 


3.067 


3.232 


6 


20 


3.482 


3.641 


4a 


5 


2.631 


2.886 


5a 


10 


3.014 


3.270 


6a 


20 


3.425 


3.688 



The possible range of DM masses for the two galaxies spans 
about half an order of magnitude. However, given the large ob- 
servational uncertainties it is the best we can do. The masses 
themselves are rather large and taken at face value would imply 
that the two dwarfs are among the most DM dominated objects 
in the observed universe. 

3.3. Comparison with observationally obtained data 

We now put our results into context and compare them with ob- 
servations. Simon & Geha (2007) measured the velocity disper- 
sion of Leo IV and derived a total dynamical mass, within their 
optical radius (97 pc), of 1.4 ± 1.5 x 10'' Mq. This implies a 
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case04a S^^^ = 26.7M^^^/pc'^ Zkpc Reso]. t = 2.37Gyr 




X [kpc] 

Fig. 3. Results of the full N-body simulation of case 4a (the case with the largest discrepancy between the two-body and full N-body 
result). In the left panel the surface-density contours of the Leo V halo are shown. The (red) cross marks the position of the centre 
of Leo IV. The inner 5 (black) contours are between 26-10 Mq pc"^, and the next 5 (blue) are between 10 and 1 and the outermost 5 
(green) contours are between 1 and 0.1. Middle panel: The solid (black) line shows the distance between the two haloes. The (red) 
horizontal line is the maximum distance from the two-body result. The right panel shows the Lagrangian radii (10, 20, 30,. ..90 per 
cent of the mass) of the two haloes; (black) solid lines for Leo IV and dashed (blue) lines for Leo V. The horizontal (red) solid line 
shows again the distance criterion from the simple code. The (red) cross shows the maximum distance in the full simulation and 
vertical dashed (red) line marks the time of maximum separation. 



M/L-ratio of 151. This value is quite similar to most of the other 
known dSph galaxies of the MW. We, therefore, use our results 
and compute the mass of our Leo IV haloes within the same 
97 pc. The resulting masses and derived M/L-ratios are shown 
in Tab.|4] Our M/L -ratios are in the range of 35-328 and encom- 
pass the results of 'Simon & Gehal (12007!). Furthermore, should 
the results of Simon & Geha (2007), which are based on very 
few stars, prove to be correct, our results mean that we can rule 
out DM haloes with high concentrations (i.e. c = 20). 

IWalker et al.l (l2009ah report a central velocity dispersion of 
cr = 2.4;^j-^ km s"' based on five stars for Leo V. Using this value 
they calculate a dynamical mass, within an adopted rh = 67.4 pc, 
of 3.3;|;25 ^ ^o- Calculating the mass within this radius in 
our models gives a range from 2.7 x 10^ Mq to 2.5 x 10^ Mq, 
again encompassing the results derived from observations by 
fWalker et al.i (i2009al) . 

Another way to compare our results with observations is by 
computing th e total mass within a 'standard' radius of 300 pc, 
as adopted bv lWalker et all ( 12009 b'). We can then infer the M/L- 
ratio inside this radius and compare o ur results with th e obser- 
vationally derived results reported by IWilkinson et al.l (l2006h . 
Given the low luminous masses of the two dwarfs, our results 
point to M/L-ratios in the order of logjQ(M/L) = 2.5-3.7 (see 
also Tab. |5]l. Despite being quite high, plotting these values 
together with observationally derived values of other dwarfs 
(Fig.|2|i we see that they follow the general trend of higher M/L- 
ratios with lower luminous masses. In fact, if we fit a line through 
the observational results, our values would intersect that linear 
fit. 

This means, not only do our results follow and confirm the 
observed trend of the known MW dSphs, we can further con- 
clude that the two satellites do not need unreasonable amounts 
of DM to form a bound pair 

3.4. Comparison to N-body simulations 

Since we now have the results of all our cases, we have to make 
sure that they still hold if we resimulate them with a full N-body 



code. Of course the full N-body simulations will differ signifi- 
cantly from the ones above, we simply want to know if our con- 
clusions remain valid. 

In the full N-body simulations, the two live haloes are in- 
teracting with each other They experience dynamical friction 
which shrinks their orbits around each other until they finally 
merge. This cannot be reproduced by the simple code but we can 
determine if the next maximum separation of the orbit is smaller 
than the extent of the other halo, as our distance criterion in the 
simple simulations requires. 

What we find is that the results differ by a few per cent (max. 
6.5%) in the simulations using a mass ratio of 1.8 between the 
haloes. The haloes get slowed down and, therefore, turn around 
at a smaller separation. 

If the haloes have the same mass and we only adopt radial ve- 
locities, the orbit overshoots the maximum distance by approx- 
imately 20 per cent. If we add the additional tangential veloci- 
ties, our restricted results differ by about 28^5 per cent, in the 
sense that the maximum separation is larger in the 'full' simu- 
lation than in the restricted one. This may seem odd given that 
dynamical friction should act in the opposite direction, however, 
there are other mechanisms at work. We see an expansion of the 
haloes as orbital energy is transformed into internal energy, fur- 
thermore, we see that the two haloes get deformed - particles 
from one halo get dragged along by the gravitational force of the 
other We give the ratio of the maximum separations between 
the full and restricted simulations in the second to last column 
of Tab. [3] (labeled 'ratio'). 

We plot, in Fig. [3] the results of case 4a, the case with the 
largest discrepancy in separation between the two haloes com- 
pared with the restricted prediction. In the left panel we see the 
contours of the Leo V halo at the time of maximum separation, 
with the cross marking the position of the centre of Leo IV. We 
see that the contours are slightly elongated towards the other 
halo and that they show a clear deformation. This deformation 
is caused by the gravitational pull of the other halo, which has 
dragged particles of the dwarf towards it. 



5 



M. Blana etal.:Leo IV & V 



In the middle panel we see the large discrepancy between the 
distance criterion (horizontal line) and the actual first maximum 
separation of the orbit. But as the total mass of a NFW profile 
only increases with the logarithm of the radius, even a large dis- 
crepancy in radius as in our case 4a amounts only to a few per 
cent error in the mass of the halo. Given the fact that our results 
span almost an order of magnitude, and the large uncertainties 
from the observations (luminous mass, distance, etc.), we claim 
that the results of the restricted code are verified. 

Finally, the right panel of Fig. [3]shows the Lagrangian radii 
of Leo IV (solid lines) and Leo V (dashed lines). We see that 
the interaction of the two haloes causes the Lagrangian radii to 
expand. At the time of the maximum separation (marked with a 
cross) the halo of Leo IV is just outside the 90 % mass-radius 
and, as shown in the left panel, is still within the expanded and 
deformed halo of Leo V. In some sense, this matches the origi- 
nal distance criterion, which said that neither halo-centre should 
leave the other halo. 



4. Discussion & Conclusions 

We have presented possible scenarios for a twin system consist- 
ing of the faint dwarf spheroidal galaxies Leo IV and Leo V. 
The simulations were performed using a simple two-body code 
to rapidly find the solutions in the vast parameter space and the 
results were verified using a full N-body code. From this we find 
the minimum DM masses required for the two galaxies to form 
a tightly bound pair. 

The parameter space is restricted by assuming two inde- 
pendent DM haloes orbiting each other. Two perfectly shaped 
haloes would only be seen before the first close passage. This is 
a strong simplification of tlie real geometry of the problem. 
But as our results show (i.e. the comparison with the real N- 
body simulations) the resulting error of this simplification is 
in the order of 5-20% and therefore much smaller than the 
mass-range of our results, stemming from e.g. the unknown 
tangential velocity. 

A further restriction is the maximum distance criterion we 
adopt. We find this criterion sensible given the satellites' large 
distances from the MW. Smaller maximum separations would 
lead to higher required masses for the system to be bound. Larger 
separations would lead to lower masses, but since the expected 
tidal radius of the system (with respect to the gravitational force 
of the MW) is of the order of our distance criterion, we feel 
confident with our choice. 

As our distance criterion is of the order of the tidal radius, we 
are able to simplify even further and treat the system of the two 
dwarfs as isolated (i.e. we do not simulate the potential of the 
MW). As our aim is to determine whether the two galaxies are 
bound now (and make no predictions about their future or past), 
we do not need to take their orbit around the MW into account. 

Regarding the relative velocity we adopt two cases. In one 
the restriction is that the measured difference in radial velocity 
is the only relative velocity the dwarfs have. In the other case the 
two satellites are given an additional tangential velocity of the 
same magnitude as the radial velocity. 

We also adopt two mass ratios. First, mass-follows-light, i.e. 
the two haloes have a mass ratio of L8 like the lumino us com- 
ponents . Second, mimicking the fact (claimed by Walk er et alj 
/2009b')) that almost all dSph galaxies reside in DM haloes of 
the same minimum mass (i.e a minimum halo mass to carry 
a luminous component), a mass ratio of 1 .0. Moreover, to span 
the full range of proposed concentrations for dwarf galaxy dark 



haloes we take three values for the concentration into account 
c = 5, 10,20. 

If we assume that the bound system consists of two DM 
haloes orbiting each other, we infer masses of about L7 - 2.1 x 
10'" M0 for the whole system. If we add an additional tangen- 
tial velocity, which cannot be observationally verified, we obtain 
X 4.2-5.4X IO'^Mq. These are indeed very high masses for the 
two faint satellites and would put them amongst the most DM 
dominated objects known. Still, these results do not infer that 
the scenario is impossible. 

Another point to take away is that if we add the same amount 
of relative velocity tangentially, i.e. increasing the total relative 
velocity by a factor V(2), the required mass more than doubles. 
This shows quite a strong dependence on the relative velocity. 
Still, if we double the tangential velocity the mass would vary 
within an order of magnitude, an uncertainty we find in our re- 
sults anyway. 

We compute the DM mass within the adopted optical radius 
of lSimon & Gehal (l2007h for Leo IV (i.e. 97 pc) and find that the 
M/L -ratios we obtain span the observationally (measured veloc- 
ity dispersion) derived results. Our results are also in agreement 
with th e measured velocity dispersion of Leo V (IWalker et al.L 
l2009al) and the inferred dynamical mass. Taking the observa- 
tions at face value, our results could, therefore, restrict the pos- 
sible concentrations of the real DM haloes, once we know their 
relative tangential velocity. 

Furthermore, we checked our results agai n st the trend 
for dSph galaxies published by IWilkinson et al.l (l2006l) . They 
give the M/L-ratios within a radius of 300 pc (also seen in 
IWalker et al.L l2009b). Our simulations predict M/L-ratios, using 
the same radius, in the range logjo M/L = 2.5-3.7. These values 
are high but encompass the predictions for faint dSph galaxies, if 
we extrapolate the known values to the magnitudes of the Leos. 

Comparing the results of our two-body code with full N- 
body simulations we find differences in the maximum separa- 
tions of only a few per cent in most of the simulations. Only 
the simulations with equal mass haloes and additional tangential 
velocity have rather large discrepancies. Because the mass of an 
NFW halo increases with radius, proportional to ln(r), the uncer- 
tainty in the masses is much lower The simple integration pro- 
gramme used cannot predict any deformations of the haloes due 
to their mutual interactions. In the full simulations we see those 
deformations and, even though the initial distance criterion is not 
fulfilled anymore, the haloes still stay within the deformations of 
the other In that sense the distance criterion is still obeyed. 

Our final remark is that we wanted to search for the necessary 
total dark matter mass of the pair of satellites to ensure that they 
are bound to each other Even though the comparison between 
the simple code and the full N-body results deviate somewhat 
from our distance criterion, they do not change the conclusions 
of the simulations. A bound pair in the restricted case is still a 
bound pair in the full simulations. Just by looking at our different 
cases (i.e. radial velocity only or radial plus tangential velocity) 
our results differ by about half an order of magnitude in total 
mass. In that respect a mass uncertainty of even 20-30 per cent 
does not change the conclusions of this paper, nor would it alter 
the inferred M/L-ratios significantly. 

Summing up, assuming that the two Leos do, in fact, consist 
of a tightly bound pair, we find their inferred dark matter masses 
to be high but still within reasonable values. Therefore, it is pos- 
sible that the two galaxies form a bound pair, making them an 
ultra-faint counterpart of the Magellanic Clouds. 
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